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I. 



INTRODUCTION 



The Yukawa potential has the form [1] 



V(r) = -- 
r 




(1) 



Here f3 is the screening parameter and r\ is the potential strength. This potential is one of the 
different types of the screened Coulomb potentials which have been studied in various areas of 
physics such as atomic physics, plasma physics, solid-state and astrophysics [2]. The readers can 
find an outline about the approaches used to study this potential in Ref. [3]. The Yukawa potential 
is a basic ground to describe interaction between charged particles in colloidal suspensions [4] and 
an important potential model in Thomas-Fermi approximation to the electron gas [5]. In the view 
of the present work, it should be noted that, with its other properties, the Yukawa potential has 
bound states only for the values of the parameter /3 below a value /3 C < 1.19 (in a.u.) [6]. The 
Yukawa potential [7-15] has also received a great deal of attention in view of the methods by 
which the potential has been studied. Some of them used to solving the wave equations are the 
1/iV-expansion [16] and shifted 1/iV-expansion methods [17], studying the potential by using the 
Raygleih-Schrodinger perturbation expansion [18], a group-theoretical approach by using the Fock 
transformation [3], the variational self-consistent field molecular-orbital method [18], the J-matrix 
method [19], a new approximation scheme proposed to study the bound states of potential [20], 
studying in terms of the hypervirial theorems [21] and a numerical solution of the Schrodinger 
equation for the present potential [22]. 

In the present work, we study the problem including also the two-component spinor wave func- 
tions in terms of the hypergeometric functions within the context of pseudospin and spin symmetries 
[23-25]. The Dirac equation with vector, V(r), and scalar, S(r), potentials has pseudospin (spin) 
symmetry when the difference (the sum) of the potentials V(r) — S(r) [V(r) + S(r)] is constant, 
which means -§:[V(r) — S(r)] = (or gjr[V(r) + S(r)] = 0). It is pointed out that these symmetries 
can explain degeneracies in single-particle energy levels in nuclei or in some heavy meson-spectra 
within the context of relativistic mean-field theories [23-35]. In the relativistic domain, these sym- 
metries were used in the context of deformation and superdeformation in nuclei, magnetic moment 
interpretation and identical bands [26]. In the non-relativistic domain, performing a helicity unitary 
transformation to a single-particle Hamiltonian maps the normal state onto the pseudo-state [27]. 
Because of these investigations, the solutions of the Dirac equation having spin and pseudospin 
symmetry have received great attention for different type of potentials such as Morse potential, 
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Eckart potential, etc. [28-32]. 

In Ref. [33] the pseudospin symmetric solutions of the Dirac equation are obtained for the 
harmonic oscillator while used the tensor interaction as a potential linear in r. In Ref. [34] spin 
and pseudospin symmetric solutions are studied for the Woods-Saxon potential by taking the tensor 
interaction as a Coulomb-like potential. In this point of view, we propose to use a Yukawa-like 
potential as the tensor interaction 



which has an attractive form. Choosing this form makes it possible to find analytical solutions of 
the present problem. 

The organization of this work is as follows. In Section 2, we briefly give the Dirac equation 
with attractive scalar and repulsive vector potentials for the cases where the Dirac equation has 
pseudospin and spin symmetry, respectively. In Section 3, we present the Nikiforov-Uvarov (NU) 
method and the parameters required within the method. In Section 4, we find an analytical 
expression for the bound states and the two-component spinor wave functions of the Yukawa 
potential by using an approximation instead of the spin-orbit coupling term. We analyze the 
problem for the cases having the Dirac equation pseudospin and also spin symmetry. We give the 
numerical energy eigenvalues for the different quantum number pairs (n, n) where we choose the 
parameter and the mass values in a.u.. The last section includes our conclusions. 

II. DIRAC EQUATION AND SPIN AND PSEUDOSPIN SYMMETRY 

The Dirac equation is basically written by using linear momentum operator, P„ = ifid^ (four- 
vector) and the scalar rest mass M. As a result, two potential couplings are used in equation. One 
coupling is a gauge invariant one to the four- vector potential A^f, t) by using P^ — > P^ — gA^ 
(g is a real coupling parameter) and the other one is to the space-time scalar potential S(f,t) 
by substitution M — >■ M + S. The " four- vector" and "scalar" terms mean that classifying the 
observable according to the unitary irreducible representation of the rotation and translation groups 
in Minkowski space-time. By taking the space component of the vector potential to vanish (A = 0) 
and writing the time component of the four- vector potential as gAo = V(r,t), then we obtain the 
so-called vector and scalar potentials, V(r) and S(r), respectively. 

The free particle Dirac equation is given (h = c = 1) 




(2) 



(i>fd li -M)*(r f ,t) = 0, 



(3) 
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and taking the total wave function as \I/(r, t) = e~ lEt ip{r) for time-independent potentials, where 
E is the relativistic energy, the above equation including also a tensor interaction, U(r), with 
spherical symmetric vector and scalar potentials is written as 



[a.P + 0M- i(3a.fU{r) + (3S{r)] ip(r) = [E - V(r)] V>(r) , 



(4) 



Here a and /3 are usual 4x4 matrices. For spherical nuclei, the angular momentum J and the 
operator K = —f3(a.L + l) with eigenvalues k = ±(j + 1/2) commute with the Dirac Hamiltonian, 
where L is the orbital angular momentum. By using the radial eigenfunctions for upper and lower 
components of the Dirac eigenfunction F(r) and G(r), respectively, the wave function is written 
as [35] 

F(r)Y^{8,4>) 
iG (r)y( 2 ) 

where Y^\d, cj)) and Y^ 2 \9, <f>) are the pseudospin and spin spherical harmonics, respectively. They 
correspond to angular and spin parts of the wave function given by 



(5) 



'2 2" 

minis 

\k\ - ^ , I = k (k > 0) ; £ 



-(k + 1) (k<0), 



(6) 



Here, Y£ me (9,(f>) denotes the spherical harmonics and mi and m s are related magnetic quantum 
numbers. 

Substituting Eq. (5) into Eq. (4) and using the followings 

(a.L)Y (2 H6, ( f>) = (K-l)Y (2 \e,<f ) ), (7a) 

(^> (1) (M) = -(«- 1)Y^(9,<P) , (7b) 

(a.f)Y( 2 H6, ( f>) = -YV(6, ( f>), (7c) 

(a.f)Y^(e, ( f>) = -Y^(e, ( f>), (7d) 

give us the following coupled differential equations 



^ + ^ - U(r) } /•'(/ ) - + \[-V(, ))G(r). 
d -- + U(r) ) G(r) = [M-E + A(r)]F(r) . 



(8a) 

- ■ , ■ l ■ - , , • (8b) 

where T(r) = V(r) — S(r) and A(r) = V(r) + Using the expression G(r) in Eq. (8a) and 

inserting it into Eq. (8b), we get a second order differential equation 



dr 2 



k(k + 1) ri\, , ,2k tt , , d Wr . . 



F(r) = - 



dT(r)/dr 



\E + M -T{r)} 



F{r) , (9) 
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where e^\r) = [E + M — T(r)] [E — M — A(r)]. By similar steps, we write the following second 
order differential equation for G(r) as 



G(r) 



dA(r)/dr 



G(r), (10) 



_[M-£ + A(r)]_ 

where e( 2 )(r) = [i£ — M — A(r)] [E + M — T(r)]. The last two equations have the following forms 

f r/ 2 k-Tk -4- 1 1 2k- J 1 

i ^2 " r I +(y " W " ^M r ) + [E + M -A][E-M - A(r)] | F(r) = , (11a) 

{ ^ " ^2^+ ( y " U(r) + ±) U{r) + [E-M-A][E + M- F(r)] } G(r) = . (lib) 

if the Dirac equation has spin symmetry which means that T(r) = A (dT(r)/dr = 0) is a constant 
and pseudospin symmetry which means A(r) = A (dk{r)/dr = 0) is a constant [25-27]. 

III. NIKIFOROV UVAROV METHOD 

The Nikiforov-Uvarov method could be used to solve a second-order differential equation of the 
hypergeometric-type which can be transformed by using appropriate coordinate transformation 
into the following form 

+ a{t)f (t) ^T + am{t) = ' (12) 

where a(t) , and a(t) are polynomials, at most, second degree, and f(t) is a first-degree polynomial. 
By taking the solution as 

= m<p® > (is) 

gives Eq. (12) as a hyper geometric type equation [36] 

-^ + ^)^r + ^T (t) = ' (14) 

where ip(t) is defined by using the equation [36] 

1 di/j(t) ir{t) 



ip(t) dt a(t) ' 

and the other part of the solution in Eq. (13) is given by 



(15) 



Mt) = W)^ [an{zMt)h (16) 

where a n is a normalization constant, and p(t) is the weight function, and satisfies the following 
equation [36] 

da(t) a(t) dp(t) 

+ -pjtj^r = T(t) ■ (17) 



The function n(t) and the parameter A in the above equation are defined as 



*(*) = ±[|a(i)-f(t)]±{± 
A = * + ^(t). 



^)-f(i) 



- cr(t) + fccr(t) 



1/2 



(18) 
(19) 



In the NU method, the square root in Eq. (18) must be the square of a polynomial, so the 
parameter k can be determined. Thus, a new eigenvalue equation becomes 



A = A„ = -n|r(t)-i(n 2 -n)^a(t) 



(20) 



where prime denotes the derivative, and the derivative of the function r(t) = f(t) + 2ir(t) should 
be negative. 

IV. BOUND STATE SOLUTIONS 

1. Pseudospin Symmetric Case 



Taking the scalar and vector potentials in Eq. (4) as 

at \ r ) s —Br t/y \ i ^ v —Br 

S(r) = e p ; V(r) = -\ e ' . 



(21) 



inserting them into Eq. (lib), using Eq. (2) and taking the following expression instead of the 
spin-orbit coupling term [37, 38, 39] 



W 1 



r 2 ^ (l_ e -0r)2 



(22) 



we obtain from Eq. (lib) 

d 2 [3 2 k(k - 1) 



-Br 



-"P i _ e -Br -Pm[E-M- A] i _ + [E-M-A][E + M] \G{r) = 



(23) 



where 771 = r] s +r/ v . Using a new variable t = 1 — e 2/3r (0 < t < 1) and the following abbreviations 



-k(k - 1) + i/(2k - v - 1) , 



= i/(-2k + 2i/) - ^ (£7 - M - A) , 



/3 



= v {-v + 1) + ^ (E - M - A) + ^ (E - M - A){E + M) . 



P 2 



(24a) 
(24b) 

(24c) 
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gives us 



d 2 



t d 



+ 



[-at-a z 2 t-alt 2 ]}G(t) = V 1 



dt 2 t(l-t)dt t 2 {l-t) 2 
Comparing Eq. (25) with Eq. (12) 



(25) 



f(t) = -t, a(t) = t(l - t) , a(z) = -a|t 2 - aid - a( , 



(26) 



Substituting this into Eq. (18), we get 



7r(t) = ^±J(^ + a 2 -k)t 2 + (-^ + a 2 + k)t + ^ + a 2 . 



(27) 



The constant k can be determined by the condition that the discriminant of the expression 
under the square root has to be zero 



& = (-\ + a 2 + k) 2 -A(± + a 2 ){\ + a 2 -k) = Q. 



(28) 



The roots of k are k\^ = —a 2 , - 2a\ =F A where A = y 7 (a| + a 2 , + a 2 )(l + 4af). Substituting these 
values into Eq. (18), we obtain 7r(t) for /ci as 



7T(t) 



1-t 



+ 



+ ai + 



A 



(29) 



and for k 2 as 



7T(t) 



1-t 



+ af + 



A 



1 n 

Now we find the polynomial r(t) from ir(t) as 



(30) 



r(t) = 1 + 2^1^? - + 2^1^? + ^== | • (31) 

so its derivative is negative. We have from Eq. (19) 

A= -al-2a\-A-\-\J^-^J=, (32) 



and Eq. (20) gives us 



2A 



A re = n I 2 + y 1 + 4af + —j===\ + n(n - 1) . (33) 
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Substituting the values of the parameters given by Eq. (24) and setting A = X n , one can find 
the following 



2^k(k-1)-±(E-M-A)(E + M) + 2^(u-l)-^E-M-A){ m + ^±^) 

+ - 2k) 2 + 4i/(l + v - 2k) + 2n + 1 = . (34) 

Thus one can find energy eigenvalues for the case of pseudospin symmetry for any K-value in 
the existence of a tensor interaction. This analytical result is consistent with the ones obtained in 
Ref. [8] for the absence of the tensor interaction. Table I presents the numerical energy eigenvalues 
for different quantum number pairs (re, k) and different values for v which makes it possible to 
see the effect of the tensor interaction on the energy eigenstates. It is seen that energy values are 
negative as stated in the theory [25-27]. Our parameter values are as follow: rji = 2.5 a.u., j3 = 0.5 
a.u., A = 0. We observe that the energy eigenvalues decrease while the principal quantum number 
increases or spin-orbit quantum number k decreases. 

In order to find the eigenfunctions, we first compute the weight function from Eq. (17) 

p(t) = t 2a ^l-t) A , (35) 

and the wave function becomes 

1 rl n 

«wM~ W^tW ^ [t " +2 °' (1 -«> b (36) 

The polynomial solutions can be written in terms of the Jacobi polynomials [40] 

¥W(t) ~p( A > 2 ^(2t-l), A>-1, ai>-l. (37) 

On the other hand, the other part of the wave function is obtained from Eq. (15) as 

^(t) = - t )(i+A)/2 _ (38) 

Thus, the radial eigenfunctions for the lower component of the Dirac eigenfunction take 

G(t) ~ t ai (1 - t)( 1+A )/ 2 p(A< 2fl i) (2t - 1) , (39) 
and the other radial component is obtained from Eq. (8a) as 

w M-E + A \ Hy J [\t 2(1 -t) J n 
+ (re + A + 2oi + lJP^t 1, 2ai+1) {2t - 1) 



(40) 



- ^P { n A ' 2ai) (2t - 1) - ^(1 - t)P^ 2ai ) (2t - 1)} . 
where we have used the property of the Jacobi polynomials as ^[Pn 9 ' r \x)] = |(re + q + r + 



2. Spin Symmetric Case 



In this case, from Eqs. (2), (21) and (22), Eq. (11a) takes the form 

— Br —fir *\ 

I Z e -fa + ^ 2 ^ + M " A l i Z e -pr + {E + M-A][E-M] jF(r) = , 



(41) 



where r\2 = r\ s — r\ v . Using the same variable gives 
d 2 t d 1 



+ 1271 ^2 ["«! " " «3^] = > ( 42 ) 



where 



dt 2 t(l-t)dt t 2 (l-t) 



-a\ = -k(k - 1) + v(2k - v + 1) , (43a) 
-a 2 . = i/(-2k + 2i/) + j (E + M - A) , (43b) 



-a^ = -v(v + 1)-^(E + M-A) + ^(E + M-A)(E-M). (43c) 
Following the same procedure we obtain the energy eigenvalues for the case of spin symmetry 



2^k(k + 1) - ±(E + M - A)(E - M) + 2^u(v + 1) + 1(£ + M - A) { m - E -^) 

+ 2k ) 2 - 4^(1 - v + 2k) - 2n - 1 = , (44) 
and the radial eigenfunctions for the upper component of the Dirac eigenfunction as 

F(t) ~ i ai (1 - t)( 1+A )/ 2 2ai ) (24 - 1) , (45) 
which gives us the other radial component from Eq. (8a) as 



,ai 1 + A k/3 (3u 



G (t) ~ " \\ "± _ 2 "i) (2t - 1) (3{l - i)(- - + + P(l - t)) 

w M-E + A \ 11 v 7 f v n t 2(1 -t) Zni frit v 7; _ 

+/3(1 - i)(n + A + 2oi + l)/^' 1+2ai) (2t - 1)1 . (46) 

The numerical energy eigenvalues for different quantum number pairs (n, k) and different parameter 
values which makes it possible to see the effect of the tensor interaction on the energy eigenstates 
are showed in Table I (for 772 = 2.5 a.u., (3 = 0.5 a.u., A = 0). We see that energy values are positive 
if the Dirac equation has spin symmetry [25-27]. We observe that the energy eigenvalues increase 
while the principal quantum number increases or spin-orbit quantum number k decreases. It could 



be interesting to study the case if we take the tensor interaction as a Coulomb-like potential. Eq. 
(2) has the following form for f3r — > 

U(r) ~ V - (1 - pr + . . .) , (47) 

which gives us an attractive Coulomb potential for \u\ < and a repulsive one for \v\ > for the 
first-order approximation. In this case, we obtain the energy eigenvalue equation for the case of 
pseudospin symmetry 



2^ K (k-1)-^2 K -1-u)-^[E 2 -M 2 -A(E + M)] + 2^(M - E + A)( m + 

- 2k) 2 - 4i/(2/c - 1 - i/) + 2n + 1 = . 

(48) 

which is valid for attractive Coulomb potential while the eigenvalue equation for the case of spin 
symmetry is written as 



2^k(k + 1) - u{2k + 1 - v) - ^ [E 2 - M 2 - A(E + M)] + 2^(M + E - A)( m + M - - ) 

+ 2k) 2 - 4^(2k + - 2n - 1 = . 

(49) 

It should be noted that the terms including the tensor interaction in Eq. (11) behave like a 
centrifugal barrier if one chooses the tensor interaction as a Coulomb-like potential. So, we could 
except that the number of bound states increase because of the contributions coming from the 
tensor terms to centrifugal barrier. We summarize the numerical results in Table II for the case 
where if we take the tensor interaction as a Coulomb-like potential given in Eq. (47). 



V. CONCLUSIONS 



We have studied the approximate bound state solutions of the Dirac equation for the Yukawa 
potential for the cases where the Dirac equation has pseudospin and spin symmetry, respectively, 
in the existence of a tensor interaction having a Yukawa-like form. We have obtained the en- 
ergy eigenvalue equations and the related two-component spinor wave functions with the help of 
Nikiforov-Uvarov method. We have presented the numerical results of the energy eigenvalues for 
the cases of pseudospin and spin symmetry in Table I to see the effect of tensor interaction on 
bound state energies and seen that this contribution is to create much strongly bound states. 
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TABLE I: Energy eigenvalues in units of a.u. of the Yukawa potential for the cases of pseudospin 

and spin symmetry (M = 1). 



pseudospin symmetry 









— U 


,, n 1 
v — U.l 


v = 


n 


K 


state 


£/ < 


h, < 


£/ < 


1 


-1 


lSl/2 


1.22470 


1.22461 


1.13610 




-2 


1P3/2 


1.46034 


1.45258 


1.23104 




-3 


14/2 


1.60000 


1.58996 


1.29851 




-4 


1/7/2 


1.69025 


1.67938 


1.34844 


2 


-1 


2S1/2 


1.13610 


1.12562 


0.84186 




-2 


2|?3/2 


1.23104 


1.21723 


0.76922 




-3 


24/2 


1.29851 


1.28329 


0.69221 




-4 


2/7/2 


1.34844 


1.33250 


0.59863 



1.0 



spin symmetry 









v = 


v = 0.1 


v = 


n 


K 


state 


£ > 


E > 


E>0 





-2 


OP3/2 


0.49589 


0.57269 


0.97850 




-3 


04/2 


1.54413 


1.55360 


1.58112 




-4 


O/7/2 


1.92856 


1.92117 


1.86174 




-5 


O59/2 


2.08474 


2.07607 


2.01207 


1 


-2 


1P3/2 






0.17028 




-3 


14/2 


1.21217 


1.26027 


1.48195 




-4 


1/7/2 


1.99593 


1.99371 


1.96986 




-5 


1/9/2 


2.19280 


2.19280 


2.17118 



1.0 
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TABLE II: Energy eigenvalues in units of a.u. for the case where tensor interaction chosen as a 



Coulomb-like potential (r?i = r? 2 = 2.5, M = 5, v = 0.5, /3 = 0.5). 







pseudospin symmetry 








spin symmetry 




n 


K 


state 


£ < 


n 




state 


E > 


1 


-1 




2.86397 





-2 


OP3/2 


— 




-2 


1P3/2 


3.36173 




-3 


0^5/2 


0.86904 




-3 


1^5/2 


3.57271 




-4 


O/7/2 


2.18996 




-4 


1/7/2 


3.67004 




-5 


O59/2 


2.78225 


2 


-1 


2S1/2 


3.40594 


1 


-2 


1P3/2 






-2 


2^3/2 


3.62177 




-3 


1^5/2 






-3 


2^5/2 


3.70932 




-4 


1/7/2 


0.03410 




-4 


2/7/2 


3.74219 




-5 


1/9/2 


1.59973 
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